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Abstract 

Given a simplicial complex K with weights on its simplices and a chain on it, the Optimal Homol- 
ogous Chain Problem (OHCP) is to find a chain with minimal weight that is homologous (over Z) to 
the given chain. OHCP has been shown to be NP-complete, but if the boundary matrix of K is totally 
unimodular (TU), it becomes solvable in polynomial time when modeled as a linear program (LP). We 
Q_i define a condition on the simplicial complex called non total-unimodularity neutralized, or NTU neu- 

-^ tralized, which ensures that even when the boundary matrix is not TU, the OHCP LP must contain an 

integral optimal vertex for every input chain. This condition is a property of the simplicial complex, and 
is independent of the input chain and the weights on the simplices. This condition is strictly weaker than 
the boundary matrix being TU. More interestingly, the poly tope of the OHCP LP may not be integral 
under this condition. Still, an integral optimal vertex exists for every right-hand side, i.e., for every input 
chain. Hence a much larger class of OHCP instances can be solved in polynomial time than previously 
considered possible. 
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Topological cycles in shapes capture their important features, and are employed in many applications from 
science and engineering. A problem of particular interest in this context is the optimal homologous cycle 
problem, OHCP, where given a cycle in the shape, the goal is to compute the shortest cycle in its topological 



class (homologous). For instance, one could generate a set of cycles from a simplicial complex using the 
persistence algorithm [15] and then tighten them while staying in their respective homology classes. The 
OHCP and related problems have been widely studied in recent years both for two dimensional complexes 



[2, 3, 5, 16, 12] and for higher dimensional instances [10, 21]. The OHCP with homology defined over the 
popularly used field of Z2 was known to be NP-hard [4]. But it was shown recently that if the homology 
is defined over Z, then one could solve OHCP in polynomial time when the simplicial complex has no 
J> relative torsion [11]. The generalized decision version of the problem considering chains instead of cycles 

(also termed OHCP) was recently shown to be NP-complete [14]. Instances that fall in between these two 
extreme cases have not been studied so far. In particular, the complexity of OHCP in the presence of relative 
torsion is not known. 

The polynomial time solvability of OHCP was shown by modeling the problem as a linear program 
(LP), and showing that the constraint matrix of this LP is totally unimodular (TU) when the simplicial 
complex does not have any relative torsion [1 1]. This connection between TU matrices and polynomial time 
solvability of integer programs (IPs) by solving their associated LPs is well known, e.g., see [20, Chap. 19- 
21]. If the constraint matrix of an LP is TU, then its polyhedron is integral, i.e., all its vertices have integral 
coordinates. Two other concepts associated with integral polyhedra that are weaker than TU matrices are 
A; -balanced matrices [6, 7] and totally dual integral (TDI) systems [20, Chap. 22] [13]. But applications of 
such weaker conditions to the OHCP LP and their potential correspondences to the topology of simplicial 
complexes have not been explored so far. 
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Our Contributions: We define a characterization of the simplicial complex termed non total-unimodularity 
neutralized, or NTU neutralized for short, which guarantees that even when there is relative torsion in the 
simplicial complex, every instance of the OHCP LP has an integer optimal solution. Under this condition, 
the OHCP instance for any input chain with homology defined over Z could be solved in polynomial time 
using linear programming even when the constraint matrix of the OHCP LP is not TU. We arrive at our main 
result by studying the structure of the OHCP LP, and characterizing several properties of its basic solutions. 
Recall that the vertices of an LP correspond to its basic feasible solutions. We prove that an OHCP LP for a 
given input chain has a fractional basic solution if and only if the OHCP LP with a component elementary 
chain, i.e., a chain with a single nonzero coefficient of 1, as input has a certain fractional basic solution. 
Using this result, we show that no OHCP LP over the given simplicial complex has a unique fractional 
optimal solution if and only if every elementary chain involved in each relative torsion has a neutralizing 
chain in the complex, i.e., when the simplicial complex is NTU neutralized. 

Our result partly fills the gap between the extreme cases of the OHCP LP solving the OHCP instance 
when the complex has no relative torsion, and the OHCP being NP-complete. The condition of a complex 
being NTU neutralized is strictly weaker than requiring the boundary matrix of the simplicial complex 
to be TU, or even to be balanced. Further, this condition is a property of the simplicial complex itself, 
and is independent of the input chain as well as the choice of weights on the simplices. Hence a much 
broader class of OHCP instances can be solved in polynomial time using linear programming than previously 
considered possible. When the simplicial complex is NTU neutralized, the linear system in the dual of the 
OHCP LP is TDI. At the same time, this case does not appear to be covered by any of the currently known 
characterizations of TDI systems. In particular, the polytope of the OHCP LP may not be integral even when 
the complex is NTU neutralized. Still, an integral optimal solution exists for every integral right-hand side, 
i.e., for every input chain. 

1.1 An Example 

We illustrate the condition of a simplicial complex being NTU neutralized by describing a set of two di- 
mensional complexes related to the Mobius strip. A 2-complex having no relative torsion is equivalent to it 
having no Mobius strip [11, Thm. 5.13]. Consider the three different triangulations of a space in Figure 1. 
In the left and right complexes, we have a Mobius strip self-intersecting at one (d) and two vertices (a, d), 
respectively, resulting in relative torsion in both cases. In the middle complex, the self intersection is along 
the edge ad, hence we do not have relative torsion. Hence the boundary matrix is TU only for the middle 
complex. Still, in the right complex, the OHCP LP has an integral optimal solution for every input chain. 

For example, consider the edge ef (shown in yellow) with coefficient 1 as the input chain. Let the edge 
weights be as follows: red, yellow, and brown edges have weight 1, green and black edges have weights 
of 0.05, and blue edges have weights of 0.10 each. The red and yellow edges are the "manifold" edges in 
the potential Mobius strip in each complex. The green edges and the single black edge am are boundary 
edges. The pairs of blue edges are boundary edges in the candidate Mobius strips, but are shared by two 
triangles each in the simplicial complex. Solving the OHCP involves pushing the heavy manifold edge(s) 
onto the light boundary edges using the boundaries of triangles. In the left complex, the optimal solution 
to the OHCP LP corresponds to all the green and blue edges with coefficients ±0.5. The OHCP solution is 
indicated in purple in the middle complex (it is the same in all three complexes). In the right complex, there 
are two integral optimal solutions to the OHCP LP - one corresponds to the purple chain, and the second 
one is shown in cyan. Any convex combination of these two chains also corresponds to an optimal solution 
of the OHCP LP, including the one made of all green and blue edges with coefficients ±0.5. This behavior 
may be explained by the presence of a disk whose boundary is an odd number of red/yellow edges, e.g., 
triangle adc, which neutralizes the Mobius strip. Adding triangle akm to the left complex makes it NTU 
neutralized. For instance, adcbnmlkja is a disc whose boundary is 9 red edges in this case. 




Figure 1 : Three triangulations of a space. The right complex is NTU neutralized, the left one is not. The 
middle complex has a TU boundary matrix. 



2 Background 

We recall some relevant basic concepts and definitions from algebraic topology and optimization. Refer to 
standard books, e.g., ones by Munkres [18] and by Schrijver [20], for details. 

Given a vertex set V, a simplicial complex K = K(V) is a collection of subsets {a C V} where a' C a 
is in K if a G K. A subset a G K of cardinality g = p + 1 is called a p-simplex. If c' C c (<r' C <r), we call 
a' a face (proper face) of a, and <r a coface (proper coface) of a'. An oriented simplex a = {vq, v±, ■ ■ ■ , v p } 
or vqVi ■ ■ ■ v p is an ordered set of vertices. The simplices o~i with coefficients a« in Z can be added formally 
creating a chain c = SjajCTj. These chains form the chain group C p . The boundary d p a of a p-simplex 
o~, p > 0, is the (p — l)-chain that adds all the (p — l)-faces of a considering their orientations. This 
defines a boundary homomorphism d p : C p — > C p _i. The kernel of d p forms the p-cycle group Z p (K) and 
its image forms the (p — l)-boundary group B p -i(K). The homology group H p (K) is the quotient group 
Zp(K)/B p (K). Intuitively, a p-cycle is a collection of oriented p-simplices whose boundary is zero. It is a 
nontrivial cycle in H p , if it is not a boundary of a g-chain. 

For a finite simplicial complex K, the groups of chains C P (K), cycles Z P (K), and H p (K) are all finitely 
generated abelian groups. By the fundamental theorem of finitely generated abelian groups [18, page 24] 
any such group G can be written as a direct sum of two groups G = F © T where F = (Z © • • • © Z) and 
T = (Z/ii © • • • © Z/tfe) with ti > 1 and ti dividing ti + \. The subgroup T is called the torsion of G. If 
T = 0, we say G is torsion-free. 

For a subcomplex Lo of a simplicial complex L, the quotient group C p (L)/C p (Lo) is called the group 
of relative p-chains of L modulo Lo, denoted C P (L, Lo). The boundary operator <9 P : C P (L) 
and its restriction to Lq induce a homomorphism 

d( L < L °):C p (L,Lo)^C p _i(L,L ). 



Cp_i(L) 



Writing Z P (L, Lq) = ker P ' for relative cycles and B p _i(L, Lq) = im 9 p ' for relative boundaries, 



we obtain the relative homology group 

Hp(£, Lq) = Z p (L, L )/B P (L, Lq). 

Given the oriented simplicial complex K of dimension d, and a natural number p, 1 < p < d, the 
p-boundary matrix of i\~, denoted [<9 P ] , is a matrix containing exactly one column j for each p-simplex a in 
if, and exactly one row i for each (p — 1) -simplex r in K. If r is not a face of <r, then the entry in row i 
and column j is 0. If r is a face of a, then this entry is 1 if the orientation of r agrees with the orientation 
induced by a on r, and — 1 otherwise. 

A matrix A is totally unimodular (TU) if the determinant of each of its square submatrix is either 0, 1, 
or —1. Hence each Aij G {0, ±1} as well. The importance of TU matrices for integer programming is well 
known [20, Chapters 19-21]. In particular, it is known that the integer linear program 

min {f T x | Ax = b, x > 0, x G Z n } (1) 

for A G z mx ™ ) b G T/ 1 can always, i.e., for every f G W 1 , be solved in polynomial time by solving its linear 
programming relaxation (obtained by ignoring x G 7L n ) if and only if A is totally unimodular. This result 
was employed to show that the OHCP for the input p-chain c modeled as the following LP could be solved 
to get integer solutions under certain conditions [11, Eqn. (4)]. 
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subject to x + — x = c + [d q ] y (2) 

x + , x > 0. 

We assume the weights Wi for p-simplices are nonnegative. Replacing y with two nonnegative variable 
vectors y + and y~, we rewrite the above LP in the following form. 



(3) 



Notice that B = [d q ], and the variable vector z T = [x +T x T y +T y T ] . Recall that xf and x i cor- 
respond to the ith p-simplex, while y+ , yj capture the coefficients for the jth g-simplex. We refer to this for- 
mulation as the OHCP LP from now on. We let P denote its feasible region, and let A = [l — / — B B~\ 
be the constraint matrix of (3). It was shown that A is TU, or equivalently, P is integral if and only if B is 
TU, which happens [11, Thm. 5.2] if and only if H p (L, Lq) is torsion-free, for all pure subcomplexes Lq, L 
in K of dimensions p and q respectively, where Lq C L. Thus OHCP can be solved in polynomial time if 
the simplicial complex is free of relative torsion. 

z is a vertex of P if it is in P, but is not a convex combination of any two distinct elements of P [20, 
Chap. 8]. A basic solution of a system of linear equations is a point in a solution space of dimension d where 
a set of d linearly independent constraints are active, i.e., satisfied as equations. If a basic solution of P is 
feasible, then it is a vertex [20, Chap. 8]. 

3 Characterizations of Basic Solutions of the OHCP LP 

Notice that P is the hyperplane defined by the equality constraints, with the only bounds being the nonneg- 
ativity constraints. We use Pa to denote the hyperplane that is P without the bounds. We use z to refer to a 
general element of R 2 ( m+n \ and call Zi an x-entry if i < 2m, and a y-entry if i > 2m. 



Definition 3.1. For any entry zi of z G R 2 ( m + n ), its opposite entry is Zj +m for i < m, Zj_ m for m < i < 
2m, Zi+n for 2m < i < 2m + n, and Zi_„ for 2m + n < i <2{m + n). We denote the opposite entry of 
Zi as z-i. Any pair of opposite entries are coefficients for the same simplex. Hence for a pair of opposite 
entries zi, z-% of z, if at least one of the two is 0, then z is concise in the zth entry, z is concise if it is concise 
in each entry. 

The following definition translates coordinates of an OHCP LP solution to the 2th row or jth column of 
[dq] , and the p- and (/-simplices the row and column represent, respectively. 

Definition 3.2. For a solution z of an OHCP LP, for any % < m, the ith p coefficient is zi — Z-i, and for any 

j : 2m < j < 2m + n, the (j — 2m) th q coefficient is Zj — z~j. 

In figures of simplices representing solutions to OHCP LPs, we generally show the p- and (/-coefficients 
of simplices, and assume that all solutions illustrated are concise. When we call a set of solutions equivalent 
we mean each has the same p- and (/-coefficients. 

For any OHCP LP, there is the unique feasible concise solution where all the y-coordinates are 0. We 
call this solution the identity solution, and denote it z 1 . For a given simplicial complex K and the constraint 
matrix A associated with its OHCP LP instances, we use z K to refer to an element of Kev(A), the kernel 
of A. For any integral z K , the set of p-coefficients of z K represent a p-chain that is null-homologous in K. 
We list some rather straightforward results from linear algebra. 

1. Any z G Pa may be written as z 1 + z K . 

2. Given z £ P A ,z = z° + z K , then z° G P A if and only if z K G Ker(A). 

3. Because A is rational, for any z K G Ker(^4), there is some scalar a > such that az K is integral. 
The following theorem is foundational to many of our later results. 

Theorem 3.3. Let z G Pa- z is a basic solution if and only if Vz G (Ker(yl) \ {0}) , 3z : z% = 0, zf~ 7^ 0. 

Proof. We prove both directions by contrapositive. Assume for some z K / 0, there is no such i. Because 
z K / 0, z + z K / z - z K . If z G Pa and z K G Ker(A), z + az x G Pa Va G K. Therefore the line 
segment L C M 2 ( m+n > defined by the two distinct end points z ± z K is contained in Pa- Hence all the 
equality constraints are active at all points in L. Consider an arbitrary inequality constraint j := Zj > 0. If 
j is active at z, then since there is no such i, zf = 0, and so j must be active for all points in L. Therefore 
any constraint active at z is active at any point in L. Therefore z cannot be a basic solution. 

Now assume z is not a basic solution. Therefore there exists some line segment L with z in its interior 
where all constraints active at z are active at all points in L. Since z G Pa, all equality constraints are 
active in L, so L C Pa- Therefore, for any other interior point z° of L we have z° / z and hence 
z — z° = z K G (Ker(A) \ {0}). Since all inequality constraints active at z are active at z°, we have 
Zi = => zf = 0, and hence zf = 0. □ 

Figure 2 illustrates nonbasic and basic solutions of the OHCP LP in a 2-complex. Orientations of 
simplices in K, coefficients of the input chains, and the p- and (/-coefficients of the solutions are shown. Note 
that the p-coefficients are the same for the second nonbasic solution and the last basic solution. Whether or 
not a solution is basic can depend on the (/-coefficients and the input chain. 



Consider the 2(m + n) x (m + 2n) matrix N 
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(entries not specified are zero). The 



columns of A^ form a basis of Ker(^4). Analyzing its structure, together with Theorem 3.3, yields the 
following results, whose proofs we omit. 



An input chain, and a basic solution. 









The same input chain, and a nonbasic solution. 

-1 





A different input chain, and a basic solution 
1 ^ ^ -1 







Figure 2: Simple examples of basic and nonbasic solutions. 

Lemma 3.4. Any basic solution of an OHCP LP is concise. 

Lemma 3.5. Any z K S Ker(A) is equivalent to a linear combination of the last n columns of N. 

Corollary 3.6. If z G (Ker(^4) \ {0}) is concise, then at least one y-coordinate in z K is nonzero. 

Lemma 3.7. Let z G Pa be a basic solution. Let z° £ Pa with z° concise in all x-entries. Let z = z° + z , 
with z K being concise in all x-entries, and for each y-coordinate j, z^j 7^ =^> Zj 7^ 0. Then z K ^ if 
and only if there exists a p-coefficient that is in z, but nonzero in z°. 

Proof. If there is a ^-coefficient that is in z, but nonzero in z°, then z K 7^ simply because z = z° + z K . 
Now assume z K / 0. Then z 7^ z°. Because z is basic, it is the only point in Pa where all entries that are 
at z are 0. Since for each y-coordinate j, z°- 7^ =^ Zj 7^ 0, there must be some x entries that are zero 
at z, but nonzero at z°. Let i be one such entry. Since i is at z but nonzero at z°, it must be nonzero in 



,K 



z . Since z and z are concise in all x-entries, — i must be in both z and z , and therefore also at z 



A' 



Therefore the p-coefficient corresponding to i and —i must be in z but nonzero in z°. 



□ 



For p = 1, Lemma 3.7 is saying that if we attempt to get to a basic solution by adding a set of triangles 
to the input chain, then adding that set of triangles must completely cancel at least one edge. Referring to 
the two basic solutions shown in Figure 2, the edge of the input chain canceled is edge ab in both cases. 



Definition 3.8. A set of vectors is linearly concise if any linear combination of the set is concise. 



The next result establishes a method for decomposing a nonbasic solution into a basic solution and a 
remainder element of Ker(A), which we will use in later analysis. 

Theorem 3.9. Let z° G Pa be a basic solution. Let z K G Ker(^4) with {zP,z, K } linearly concise. Let 
z = z° + z . Then z is a basic solution if and only if there do not exist z , z satisfying the following 
properties: 

1. z c + Z D = Z K . 

2. z c ,z D e Ker(A). 
J. z D / 0. 

4. {z°, z , z } w linearly concise. 

5. z° + z c = z 1 w a oas/c solution. 

6. For each y-coordinate j, zS ^ =^> 2j 7^ 0. 

7. For eac/i x-coordinate i, zf 7^ =^> Zj 7^ 0. 

Proof. Suppose there exists such a decomposition of z K . By Property 2, and that z° G Pa, we have that 
z G Pa- Therefore if z is a basic solution, then the vectors z, z 1 , z D satisfy the conditions for z, z°, z^ in 
Lemma 3.7. Because of Property 4, {z°, z^, z c , z 15 , z 1 , z} is linearly concise. Therefore (Vx-coordinate i , 
zi = z-i =^> z\ = z}_j) <^=^ (Vx-coordinate i , z% = =^> z\ = 0) -4=> (Vx-coordinate i, Zi = 
=> zf = 0) -4=>- (Vx-coordinate 2 , zf / =^> z% / 0). So by Lemma 3.7, z is not a basic solution. 
Now suppose z is not a basic solution. {z°, z^", z} is still linearly concise. Construct z D and find z 1 
using the following algorithm. 

1. Let z D = 0, z 1 = z. Then {z°, z^", z 15 , z 1 } is linearly concise. 

2. z 1 must be in P4, and is not a basic solution. By Theorem 3.3, Bz" 6 (Ker(^4) \ {0}) where 
zf / =^> z\ / 0. Because {z°, z K , z 15 , z 1 } is linearly concise, {z°, z^, z D , z 1 , z w } is linearly 
concise. 

3. By Corollary 3.6, zf / for some y-coordinate i. Find i such that zj^ / 0, j > 2m =>■ \z\jzf\ < 
\z 1 /z N \ 

4. Let a = z}/zf. 

5. Let z D = z D + az^, z 1 = z 1 — az^. Because we may add any linear combination of a set of linearly 
concise vectors to the linearly concise set, {z°, z K , z D , z 1 , z w } is still linearly concise. 

6. IF z 1 is not a basic solution THEN LOOP to Step 2. 

7. STOP. 

Because if is finite, and we make at least one y-entry zero that was nonzero in z 1 in each loop, and do not 
make any zero y-entries in z 1 nonzero. Hence by Theorem 3.3, this algorithm must eventually terminate. 
More precisely, it must terminate after at most n iterations. By our criteria of choosing z N in each loop, 
z D , z 1 , and z*- 7 = z^ — z D satisfy all the criteria of the theorem. □ 

The next Lemma 3.10 describes necessary conditions to transform a nonbasic solution to a basic one. 



Lemma 3.10. Let z° G Pa be concise. Let z = z° + z 1 where z°, / =^- z} = \/j > 2m. If z is a 

basic solution in Pa, then for each z 6 (Ker(yl) \ {az 1 }) (a G R) where zf ^ =^> zf 7^ 0, there 

z K z K 
must be two x-coordinates r and s where z r = z s = 0, zfr , zf 7^ 0, and -y ^ -Ar-. Furthermore, if O r 

z r z s 

and O s are the OHCP LPs with input chains where the only nonzero coefficients are r and s, respectively, 
with these coefficients equaling those in z°, then z 1 + z is a basic solution to O r or O s where z is the 
solution with {z 1 , z 1 '} linearly concise and equivalent to the identity solution for O r or O s , respectively. 

Proof. By Theorem 3.3, the existence of z K G (Ker(A) \ {az 1 }) where zf / =^> zf / 0, implies z° 
is not a basic solution. Assume there is such a z K with no such r and s. By Zj / =^> z} = Vj > 2m, 
we get zf / =^> Zj / \/j > 2m. So if there is no x-coordinate r where zf / 0, z r = 0, then 
z cannot be a basic solution. Let ^ be the set of x-coordinates where zf / 0, z r = 0, and assume & is 
nonempty. Because there is no such r and s, there is some a G M such that {zf / z}) = a\/r G ^. If 
z G Pa, then z 1 G Ker(A). So then z K - az 1 is in Ker(^l), and Zi = =^> zf" - a^ 1 = Vz < 2m. 
Because z K / az 1 , z^ — az 1 / 0. Therefore by Theorem 3.3, z is not a basic solution. 

Let z 1 ' be the solution with {z 1 , z 1 '} linearly concise and equivalent to the identity solution for O r . If 
z 1 ' + z 1 is not a basic solution to O r , then we may decompose z 1 into z c + z D according to Theorem 3.9. 
If z D does not bring z® to zero in our original OHCP LP O, then z cannot be a basic solution of O because 
all nonzero coefficients of z D will be nonzero in z. Even if it does bring z® to zero, z° + z c cannot be a 
basic solution because of the first part of this lemma, and because z D does not bring any coefficients of z c 
to zero, by the first part of this Lemma, z° + z c + z D cannot be a basic solution unless there is another r and 
s satisfying all qualities of the lemma. A symmetric argument holds replacing r with s and vice versa. □ 

We illustrate the Lemma in Figure 3. In this example, r and s correspond to the edges be and ea, 
respectively, of z 1 . Note that both of these edges have coefficients of —1 in z 1 , but one has a coefficient of 
1 in z K . Therefore the inequality of the ratios specified in Lemma 3.10 holds; one ratio is 1, and the other 
—1. Isolating the edges corresponding to r and s as inputs with coefficients taken from z° will show that the 
other requirements of Lemma 3.10 are also met. 

The remaining results of this section describe the relationship between the existence of (non)integral 
basic solutions of an OHCP LP, and the existence of (non)integral vertices. 

Lemma 3.11. Let z° be concise. For any z, there exists a z' such that {z', z } is linearly concise, and z and 
z' are equivalent. 

Proof. For each i such that z^ / 0, z°^ i / 0, subtract Zi from both zi and z_j. The result will be both 
equivalent to z', and form a linearly concise set with z°. □ 

Lemma 3.12. z is a basic solution if and only ifz is concise, and each z' that is concise and equivalent to 
z is a basic solution. 

Proof. Assume z is a basic solution. Then by Lemma 3.4, it is concise. For any z' that is also concise and 
equivalent to z, we transform z to z' by taking each Zi where zi / z\, and subtracting z% from both zi and 
z-i. For each pair of subtractions, we make exactly one inequality constraint active, namely zi > 0, and 
exactly one inactive, z_j > 0. For any linearly independent set of constraints containing Zi > 0, if we 
replace this constraint with z-i > 0, the result must again be a linearly independent set. 

If each z' that is concise and equivalent to z is a basic solution, and z is concise, then a similar logic 
holds to show z is a basic solution. □ 

Corollary 3.13. Ifz is a basic solution, then there is a unique vertex z' that is equivalent to z. 




z = z° + z 1 = (z 7 + Z K ) + z 1 




@ (7> 




Figure 3: A Mobius strip illustrating Lemma 3.10. 



Proof. For each i with Zi < 0, subtract Zi from both Zj and z_j. The result will be both concise and 
equivalent to z, and so by Lemma 3. 12 will be a basic solution. It will also be nonnegative, and therefore a 
vertex. Also for each i, there is one value we may add to both z% and z_ j such that the result will be a vertex 
and equivalent to z. □ 

Corollary 3.14. Let z be concise. Then z z's integral if and only if each z' ?/to w concise and equivalent to 
z jj integral. 

Proof. We transform z to z' using the same method as in Lemma 3.12. The result holds by the closure of 
integers under addition. □ 



4 Fractional Solutions to the OHCP LP, and Elementary Chains 

Consider the special case of OHCP where the input chain c is the elementary p-chain ej for some i < m, 
i.e., the zth p-simplex has a coefficient of 1 in the chain while all other entries are zero. We refer to this 
instance of the OHCP LP as OHCPj, and its feasible region as Pj. We analyze the relationship between the 
existence of nonintegral values in basic solutions of the OHCP LP, and the existence of nonintegral values 
in a basic solution of OHCPj for some i. 

Lemma 4.1. Let a 6 (R \ {0}). z is a basic solution of the OHCP LP O with input chain c if and only if 
az is a basic solution of the OHCP LP O' with input chain cue. 



Proof. If a 7^ 0, then multiplying z by a does not change which entries are nonzero. Therefore the set of 
active inequality constraints are the same at z and az. And if we multiply both a solution and the input 
chain by the same scalar, the set of active equality constraints cannot change. Therefore the set of active 
constraints is the same for z in O and for az in O'. Therefore z is a basic solution to O if and only if az is 
a basic solution to O'. □ 

Theorem 4.2. Let z £ Pa be a basic solution to the OHCP LP, with {z, z 1 } linearly concise. There exists 
some matrix Z such that the columns of Z form a linearly concise set, each column Zi of Z is a basic 
solution in (Pi) a, an d Zc = z. 

Proof. Assume z is a basic solution in Pa- Begin constructing Z by first setting each column Z\ to the 
vector that is linearly concise with z 1 , and equivalent to the identity solution of OHCPj. Then Zc = z 1 . By 
Theorem 3.3 and Lemma 3.5, the identity solution to any OHCP LP is a basic solution. So by Lemma 3.12, 
each Z{ is a basic solution to OHCPj in (Pi) a- Also, by Lemma 4.1, aZi is a basic solution to the OHCP 
LP with input chain Cjej. 

Let z^ = z — z 1 . If z K / 0, distribute it among the columns of Z according to the following algorithm. 

1. Letz° = z J . 

2. Since {z, z } is linearly concise, z K is concise. By Lemma 3.7, there is some x-coordinate i such 
that zf / 0, Zi = 0. It must also be true that cy / where i' is either i or —i. Note that Zy is equal 

to[e£ T T T ] T 

3. IF cyZy + z K is not a solution to the OHCP LP with input chain cyey, THEN 

(a) Construct z c and z D according to the algorithm in Theorem 3.9. Note that z° and z K represent 
the same vectors in Theorem 3.9 as they do here. However, the difference between z of this 
Lemma and z of Theorem 3.9 is z 1 — cyZy. 

(b) z c cannot be 0. Otherwise, z D would be z K , and so by our choice of i in step 2, would not 
satisfy Property 7 of Theorem 3.9. Also recall from Theorem 3.9 that all vectors referred to in 
this algorithm form a linearly concise set. Add —z c to Zy. 

(c) Set z K := z D . By Property 3 of Theorem 3.9, we still have z K / 0. 

(d) Set z° := z° + z c . {z,z } is still linearly concise. And because z c satisfies Property 6 of 
Theorem 3.9, we still have that for each y-coordinate j, z°- / =^> Zj / 0. Also, by 
Property 7, for any x-coordinate i where zf / and Zi = 0, cy ^ where i' is either i or — i. 
This is because the only x-coefficients that are not zero in the z of Theorem 3.9, but are zero in 
the z of this theorem must be nonzero in z 1 . 

(e) LOOP to Step 2. Note that the next i chosen in Step 2 cannot be the same as any previous is 
chosen, since by Lemma 3.7, and our previous choices of i, z® must be for each previous i. 

4. Add ^-z K to Zy . 

5. STOP. 

Because each z c is nonzero, and K is finite, this algorithm must terminate, giving us the desired Z. □ 

Consider the case of p = 1. Because a basic solution z cannot be decomposed as in Theorem 3.9 the set 
of triangles of z form a union of 2D spaces each of which must be connected to at least one edge r of the 
input chain c. Then each of these 2D spaces must be a basic solution to the OHCP LP where r has the only 
nonzero coefficient in the input chain c T , and the coefficient of r is the same in c and c r . By Lemma 4.1, 
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we may scale each of these basic solutions as necessary to be basic solutions of elementary chains. These 
basic solutions of elementary chains, adjusted as necessary to form a linearly concise set, are columns of 
Z. Note that the construction of Z may not be unique, and many of the columns of Z may be equivalent to 
identity solutions. 

We may refer to Figure 2 and think of the edges of the last input chain as three different input chains 
consisting individually of ab, be, and ca. We may then decompose the basic solution shown into solutions 
equivalent to identity solutions for be and ca, and a basic solution to ab. 

Lemma 4.3. For a given complex K, there is an OHCP LP with integral input chain c that has a nonintegral 
basic solution if and only if there is some i such that the ith coefficient is nonzero in c, and OHCPi has a 
nonintegral basic solution. 

Proof. If there is some i such that OHCPi has a nonintegral basic solution, then ej is an integral input chain 
that has a nonintegral basic solution. If z is a nonintegral basic solution to the OHCP LP with integral 
input chain c, then by Theorem 4.2, z is the sum of n vectors each of the form Cj {z 1 ) , . + z c , where 
•zP / =^> Ci ^ 0. Since integers are closed under addition, one of these terms must be nonintegral. 
Since c, and (z 7 ) , . both must be integral, there must be some z c that is nonintegral. Therefore some ^z c 
is also nonintegral, and so (z 1 ) , . + — z c is nonintegral. By Theorem 4.2, each (z 1 ) , . + — z c is a basic 
solution to OHCPj. Therefore some {z 1 ) , . + ^z c is a nonintegral basic solution to OHCPj. Furthermore, 
the ith coefficient of c must be nonzero, otherwise this sum is undefined. □ 

Lemma 4.4. Let z° be concise and be in the hyperplane Pa of the OHCP O with input chain c, with z a 
basic solution to the same OHCP in Pa where for each y-coordinate j, Zj = =^> z^ = 0. Let y° be the 
2(ra + n)-vector with all y-coefficients equal to those ofz°, and all x- coefficients 0. Then z — y° is a basic 
solution of the OHCP 0° with input chain c° in {Pa) where [(c°) T ] T is equivalent in all x-coordinates 
toz°. 

Proof. Since z and z° are both in Pa, the difference between them is in Ker{A). Since [(c°) T ] T and z° 
are equivalent in all x-coordinates, z° — y° is equivalent to the identity solution for O . Therefore z° — y° 
and z — y° are both in {Pa)°- 

We show that z — y° is a basic solution of O by contrapositive. Let z* represent z — y°. Suppose z* is 
not a basic solution of 0°. Then by Theorem 3.3, 3z K G (Ker(yl) \ {0}) where z* = =^ zf = 0. z 
and z* agree in all x-coordinates. We also have for each y-coordinate j, Zj = =^ z, = =^> y' = 
=^ z* = =^ z^- = 0. Therefore we have Zi = =^ zf = for all i. Therefore by Theorem 3.3, 
z is not a basic solution of O. □ 

For p = 1, due to Corollary 3.13, Lemma 4.4 is saying that if we have a set of edges x that is a vertex 
for some OHCP LP with input chain c, and we transform c to x by adding triangles one at a time, or at least 
not eliminating triangles previously added, then x will be a vertex to any OHCP LP that has as input any of 
these intermediate edge sets we get at each step of this transformation. 

Lemma 4.5. For a given complex K, there is an OHCP LP with integral input chain c that has a nonintegral 
basic solution if and only if there is some i such that OHCPi has a basic solution where all y-coordinates 
that are nonzero are nonintegral. 

Proof. Let z be a basic solution with nonintegral coefficients for an OHCP LP with integral input chain c. 
By Theorem 4.2, there is some matrix Z such that the columns of Z form a linearly concise set, each column 
Zi of Z is a basic solution in {Pi) a, and Zc = z. By a similar logic as Lemma 4.3, one of these columns is 
nonintegral, and this column is of the form {z 1 ) , . + ^z c ' . For this z c ', let ^ be the set of y-coordinates 
with nonzero integral coefficients in — z c . 
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If J* is empty, then we have our desired result. If fl is nonempty, decompose z c into its linear combi- 
nation of basis vectors of Ker(^4) that are equivalent to the last n columns of the matrix N of Lemma 3.5. 
Let z° be the sum of e±i and the components of this linear combination with nonzero element in ^ . Note 
that z° is integral. Let c° be the input chain where [(c°) T ] T is equivalent in all x-coordinates to z°. Let 
y° be the 2(m + n)-vector such that j G J! =^> \p- = (1 / a) zf , j ^ J? =^> y^ = 0. Then by 
Lemma 4.4, (z 7 ) , . + (1/cj)z c — y° is a basic solution to the OHCP LP O with input chain c°, and all of 
its nonzero y-coefficients are nonintegral. 

Now let z = (z 7 ) , . + (1/cj)z — y°, let c = c°, and apply the same logic as above. Note that with each 
iteration of this process, we lessen the number of nonzero y-coefficients; the set of nonzero y-coefficients of 
z contains the set of nonzero y-coefficients in (1/cj)z , and this set decreases in size each round. Because 
K is finite, this process must eventually terminate. □ 

Remark 4.6. Lemmas 4.3 and 4.5 allow us to concentrate on the easier to analyze case of elementary input 
chains for the OHCP LP in order to arrive at our main results. 

5 Projections Onto the Space of j>Simplex Coefficients 

Call the space of x-variables X (= R 2m ) . We study the projections of the OHCP LP, P, P A , Ker(A), basic 
solutions, and vertices onto SC . For any of these objects Q, let Q\ %■ represent the projection of 17 onto SC . 
Since the y-variables do not appear in the OHCP LP objective function, a vertex of P\% must be optimal. 
We extend the definition of concise to (2m)-vectors in the natural way. We also rephrase the definition of a 
basic solution below. 

Definition 5.1. Let <?f be the set of constraints of O that are not orthogonal to SC , i.e., the set of constraints 
with a nonzero x-coefficient. For any x G X , let <^ x be the set of elements of *£ active at x. Then x is a 
basic solution ofO\% if and only if there is no other point in X where all elements of ^ x are active. 

Note that if x fails Definition 5.1, then <jf x must be active for an entire affine space of some dimension 
d > 1. To say x € X is feasible in 0\x is equivalent to saying x G P\%- We justify Definition 5.1 and 
how it relates to vertices of P\ sc with the following lemma. 

Lemma 5.2. x G X is a basic feasible solution of ' 0\% if and only if x is feasible, and not a convex 
combination of any two distinct elements of P\sc- 

Proof. We prove both directions by contrapositive. If x is infeasible, then clearly it is not a basic feasible 
solution. Suppose x is a convex combination of two distinct elements x 1 and x 2 of P\%. Then x = 
Ax 1 + (1 — A)x 2 for some A € (0, 1). Since x x ,x 2 G P\$r, they are nonnegative, and so the set of 
nonzero coefficients of x is the union of the sets nonzero coefficients of x 1 and x 2 . Therefore any inequality 
constraints in ^ active at x must also be active at x 1 and x 2 . Since P is convex, P\x is convex. Since 
x : ,x 2 G P|^,x G P\gc- Therefore all equality constraints are active at all three points. Therefore x does 
not satisfy Definition 5.1 as a basic solution of 0\x- 

Now suppose x is feasible, but does not satisfy the conditions specified in Definition 5.1. Let x' be 
another point in SC where all elements of <jf x are active. All equality constraints are active at x, and so are 
also active at x'. Therefore x' G (P\%-) A . Let L be the line in X containing x and x'. Then L C (P\%r) A . 
Any point in L may be expressed as ax.' + (1 — a)x for some a G R. Choosing a value for a defines a 
point. Because P\<% is convex, there is at most two values for a such that ax' + (1 — a)x is on the boundary 
of P\x- Since x is defined by a = 0, and x is feasible, at most one such value is positive, and at most one 
such value is negative. 
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If there is no such value for a, then L C P\ %. Then choose values 1 and —1 for a, to define x 1 and x 2 . 
If there is only one such value a, then choose a and —a to define x 1 and x 2 . If there are two such values 
«i, CK2, then choose min{|ai|, \a2\} and — (min{|ai|, |«2|}) to define x 1 and x 2 . In any case, both x 1 and 
x 2 areinP|, r ,andx= (l/2)x 1 + (l/2)x 2 . D 

Lemma 5.2 shows we may define vertices of P\ % the same way as in P. We now show results for basic 
solutions of 0\f£ and vertices of P\% that parallel many of our previous results. Proofs are omitted where 
the logic is a natural parallel of these previous results. 

Corollary 5.3. Letx € (Pa) \x- ~x~ is a basic solution of 0\ %■ if and only if\/x. K G (Kei(A)\_%- \ {0}) , 3i : 
Xi = 0,x? ^0. 

Lemma 5.4. Any basic solution xofO\%- is concise. 

A significant difference in basic solutions of the projection from basic solutions in the whole solution 
space is that in the projection, because they are determined only by p-coefficients, the choice of input chain 
on a complex has no impact on whether or not a solution in (Pa) \x is basic in the projection. Referring 
back to Figure 2, the last solution is not basic in the projection. 

Lemma 5.5 formalizes this idea of the input chain not mattering, with some added context useful for our 
main result. It is a parallel to Lemma 4.4, but becomes simpler in the projection space. 

Lemma 5.5. Let z° be concise and be in the hyperplane Pa of the OHCP O with input chain c, with z a 
basic solution to the same OHCP in Pa and z\% a basic solution to 0\x- Then z\% is a basic solution of 
the OHCP O \$r with input chain c° in (Pa)°\x where [(c°) T ] T is equivalent to z°\$r. 

Proof. Since z an z° are both in Pa, the difference between them is in Ker(yl). Since [(c°) T ] T and z° 
are equivalent in all x-coordinates, z°\%- is equivalent to the identity solution for 0°. Therefore z°|^ and 
z\% are both in (Pa)°\%- 

We show that z\%- is a basic solution of 0°|x by contrapositive. Let x = z| ^. Suppose xis not a basic 
solution of 0°\%-. Then by Corollary 5.3, 3x K £ (Ker(A)\%- \ {0}) where x t = => xf = 0. But 
then by Corollary 5.3, z\% is not a basic solution of 0\ % . □ 

Lemma 5.6. For any basic solution x G (Pa) \x °f 0\^, there is a basic solution z € Pa of O where 
z\x = x. 

Proof. Since x e (Pa) \sc, there is some z' where z'|x = x, and z' G Pa- Since z' and x agree in all 
x- variables and x is a basic solution, for any z K where Zi = =^ zf = 0, z K \% must be 0. If such 
a z K exists, we may use the algorithm of Theorem 3.9 to arrive at a basic solution z of O in Pa where 

z| x = x. □ 

Corollary 5.7. If for a given complex K, all vertices of any OHCP LP are integral, then all vertices of any 
projection of an OHCP LP onto & must be integral. 

Proof. We prove by contrapositive. If there exists some 0\x with a nonintegral vertex x, then by Lemma 5.6 
and Corollary 3.13, there is a vertex z of O where all x-coordinates of z and x agree. Therefore z is nonin- 
tegral. □ 

Corollary 5.8. Let x € (Pa) \x be a basic solution of 0\%. Let x° £ (Pa) \x with x° concise. Let 
x = x° + x , with x being concise. Then x ^ if and only if there exists a p-coefficient that is zero in 
x, but nonzero in x°. 
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Corollary 5.9. Let x° G {Pa) \% be a basic solution ofO\%. Let x. K G Ker(^4) | % with {x°, x^} linearly 
concise. Let x = x° + x . Then x is a basic solution if and only if there do not exist x , x satisfying the 
following properties: 

1. x c + x D = x x 

2. x c ,x D eKerL4)|<r. 

3. x D / 0. 

4. {x°, x K , x B } j's linearly concise. 

5. x° + x = x 1 w a ^as/c solution to 0\%. 

6. x? ^ o => Si / Vi. 

Proof. We prove both directions again by contrapositive. The first direction follows in the same way as 
Theorem 3.9, replacing z with x and Lemma 3.7 with Corollary 5.8. Then suppose x is not a basic solution. 
{x°, x. K , x} is still linearly concise. Construct x D and find x 1 using the following algorithm. 

1. Let x D = 0, x 1 = x. Then {x°, x^, x 15 , x 1 } is linearly concise. 

2. x 1 must be in {Pa) \x, and is not a basic solution. By Corollary 5.3, 3z^ G (Ker(^4)|x \ {0}) 
where xf / =^> x\ / 0. Because {x°, x. K , x D , x 1 } is linearly concise, {x°, x^, yp , x 1 , x^} is 
linearly concise. 

3. Find i such that x^ ^ =>■ \x) lxf\ < \x]/x^\. 

4. Let a = xj/xf . 

5. Let x-° = x D + ax^jX 1 = x 1 — ax^. Because we may add any linear combination of a set of 
linearly concise vectors to the linearly concise set, {x°, x x , x D ,x 1 ,x ff } is still linearly concise. 

6. IF x 1 is not a basic solution of 0\ x THEN LOOP to Step 2. 

7. STOP. 

Because K is finite, and we make at least one entry zero that was nonzero in x 1 in each loop, and do 
not make any zero entries in x 1 nonzero, then by Corollary 5.3, this algorithm must eventually terminate 
after at most m iterations. And as in Theorem 3.9, x D , x 1 , and x c = -x. K — -x. D satisfy all criteria of the 
Corollary. □ 

Corollary 5.10. Let x° G {Pa) \x be concise. Let x = x° + x 1 . If~xisa basic solution in {Pa) \x> 
then for each x^ G (Ker(A)| <% \ {ax 1 }) for a G M, where xf 7^ =^> xf 7^ 0, there must be two 
coordinates r and s where x r = x s = 0, xf , xf / 0, and {x r /x r .) / {x s /x s ). Furthermore, if{O r ) \% 
and {Og) \% are the projections ofOHCP LPs with input chains where the only nonzero coefficients are r 
and s respectively, with these coefficients equaling those in x°, then x 1 + x is a basic solution to {O r ) \x 
and {O s ) \% where x 7/ is the solution with {x 1 , x /; } linearly concise and equivalent to the identity solution 
for (O r ) \x and {O s ) \%, respectively. 

Corollary 5.11. x is a basic solution of 0\g; if and only ifx is concise, and each x' that is is concise and 
equivalent to x is a basic solution. 

Corollary 5.12. If x is a basic solution, then there is a unique vertex x' that is equivalent to x. 
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Corollary 5.13. Let x be concise. Then x is integral if and only if each x' that is concise and equivalent to 
x is integral. 

Theorem 5.14. For a given complex K, there is a nonintegral vertex z of some OHCP LP O with integral 
input chain c and polyhedron P where z\%- is a vertex ofP\x if and only if there is some i such that OHCPi 
has a nonintegral vertex z' such that z'\^ is a vertex ofPf\% where the ith coefficient of c is nonzero. 

Proof. One direction of the if and only if is trivially true: if there is some i and OHCPi, then the more 
general case of OHCP LP follows immediately. We prove the other direction by contrapositive. If there is 
no such % where OHCPj has a nonintegral vertex z', then by Lemma 4.3 and Corollary 3.13, there can be no 
nonintegral vertex z of any OHCP LP. 

Now suppose there is an i where OHCPj has a nonintegral vertex z', but no such i where z'\& is a 
vertex of P%\sc- Then for any nonintegral vertex z of an OHCP LP, by Theorem 4.2 and Lemma 4.3, some 
column z' of Z with Zc = z is a nonintegral basic solution of OHCPj for some i. If we construct Z using the 
algorithms of Theorem 4.2 and Theorem 3.9, then by Condition 7 of Theorem 3.9, any nonzero p-coefficient 
of z' is nonzero in z. 

Since z'|#- cannot be a vertex, it cannot be a basic solution of OHCPj \<%. So by Corollary 5.9, there is 
some yp G Ker (A) \ % where every nonzero coefficient of yp is nonzero in z' | %■> Since all these coefficients 
are ^-coefficients, all the nonzero coefficients of x D are nonzero in z| jr. So by Corollary 5.9, z\x cannot 
be a basic solution to 0\sc, and so is not a vertex of P\%. □ 

Corollary 5.15. For a given complex K, there is an OHCP LP with integral input chain c and polyhedron P 
that has a nonintegral vertex z where z\% is a vertex of P\% if and only if there is some i such that OHCPi 
has a nonintegral vertex z' with all nonzero y -coefficients nonintegral where z'\% is a vertex of Pi\%. 

Proof. The result follows from Lemmas 4.5 and 5.5, and Corollaries 3.13, 3.14, 5.12, and 5.13. □ 

6 Minimally Non Totally-Unimodular Submatrices of [d q ] 

A minimally non totally -unimodular (MNTU) matrix is a matrix M that is not totally unimodular, but every 
proper submatrix of M is totally unimodular (also referred to as almost totally unimodular matrices [1]). 
The following properties of an MNTU matrix M are already known [1, 22]. 

1 . A matrix is not totally unimodular if and only if has an MNTU submatrix M. 

2. det(M) = ±2. 

3. Every column and every row of M has an even number of nonzero entries, i.e., M is Eulerian. 

4. The sum of the entries of M is 2 mod 4. 

5. The bipartite graph representation of M is a chordless (i.e., induced) circuit [8]. 

The bipartite graph representation [1, 8] a submatrix M of [d q ] has a vertex for each row and for each column 
of M, and an undirected edge for each nonzero entry M^ connecting the vertices for row % and column j. 
Notice that each edge connects a row vertex, or p- vertex, with a column vertex, or g-vertex. A circuit C in 
a weighted graph is b-odd (b-even) if the sum of the weights of the edges in C is 2 mod 4 (0 mod 4). The 
quality of C being b-even, b-odd, or neither is called the b-parity of C. The following theorem characterizes 
this bipartite graph as a circuit. 
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Theorem 6.1. Given a circuit C that is the bipartite graph representation of an MNTU submatrix M of [d q ], 
and a set of flags placed on an arbitrary subset of the q-vertices ofC, there exists a traversal ofC such that 
each portion of the traversal ofC between two consecutive flags is induced. 

Proof. Recall that for any circuit C, if an edge h G C is a potential chord of a subgraph of C, then both 
of its end points are of degree 4 or more in C. If there is no path in C between two flags, then every path 
between them must contain both end points of a chord, and we have a set of paths like one of the graphs in 
Figure 4. The graphs shown are abstractions of C in the case where the end points of h\ and /12 all have 
degree 4, which is the minimum possible degree for these vertices. Graph C\ is the case showing the four 
half-paths, and the other graphs show the possibilities of how these half-paths can connect. 



a 3 a 






Figure 4: Abstract representations of C. 



In graph C2, though any path between the two flags contains both ends of a chord, we may still traverse 
the entire graph in such a way that each portion of the traversal between flags is induced. This can be done 
by traversing a\ and 05 immediately before or after hi, and traversing 0,3 and a,Q immediately before or after 
/12. In this way, all paths between end points of a specific chord that do not contain a flag are traversed 
consecutively. 

In graph C3, we may also traverse the graph in such a way that each portion of the graph is induced. If 
we traverse either 05 or ae as soon as possible, then we have a path between the flags that does not contain 
both end points of any chord. We may then traverse that remainder of the graph, which is also induced. 

In graph C4, we cannot traverse the graph without one of the portions between the flags not being 
induced. This graph may be decomposed into four cycles: a\ and hi, a-i and a§, 03 and /12, as well as 04 
and G15. Note that of these four pairs of paths, there must be at least one pair where neither path of the pair 
is a chord. Otherwise, C would have a cycle of four edges. And because no two distinct g-simplices share 
two distinct p-faces, this is impossible in a submatrix of [d q ]. 

Suppose without loss of generality that neither 04 nor 05 is a chord. Then these two paths form a cycle. 
If this cycle is not induced, that means it contains both end points of a potential chord not shown, but not the 
chord itself. Then we alter the two paths so that whenever we encounter an end point of a potential chord, 
not in either path, we always choose to cross the chord. Each altered path must still have the same end points 
shown in C4, or C4 would not apply. 
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This gives us a cycle Y that is induced. While 04 or 05 may contain a potential chord, it is still true that 
neither can be a chord. Therefore C \ Y is also induced. Because C is b-odd, either Y or C \ Y must be 
b-odd. But this contradicts M being minimal. Therefore graph C\ cannot occur. 

If we suppose that one or more of the end points of the chords is of degree more than 4, each of these 
points still must be of even degree. If any added paths connect diagonally, then we have the same case as 
graph C 3 . 

If no added paths connect diagonally, then divide the paths of the graph into four sets: those connecting 
the tops of two different chords, those connecting the bottoms of two different chords, and two sets con- 
necting end points of the same chord. If any added path loops back to its starting vertex, ignore it for now. 
Because all four end points must be of even degree, if any one of these sets contains an odd number of paths, 
then they all must. This is equivalent to case Ci- If all four sets contain an even number of paths, then this 
is equivalent to case C4, and cannot occur. Note that any paths we have ignored that loop back to the vertex 
where they began do not affect the parity of these other four sets of paths. 

Also note that paths shown in these abstractions may cross each other or themselves in ways not shown, 
but this will not affect our results, as explained below. 

• If two paths that both connect the same end points cross, then because we are only discussing the 
existence and parity of the number of paths between end points, our argument is unaffected. 

• If a path connecting the two top points of the chords crosses a path connecting the two bottom points, 
then Case C3 applies. 

• Case C3 also applies if a path connecting the ends of the same chord crosses a path connecting the 
ends of the other chord. 

• If a vertical path crosses a horizontal path, or if a path that starts and ends at the same point crosses 
paths from two different sets, it is possible none of the graphs shown apply. But there still must be a 
path between the flags that does not contain both end points of any chord. 

Also note that the flags may actually be placed at an end point of a chord. However, because the flags 
may only be placed at g-vertices, and C is bipartite, they cannot be placed at both ends of the same chord. 
Therefore the placement of flags will also not affect our results, and we may abstract this placement as 
shown. 

As a final step to show our result that a full traversal exists, simply cut out a hole in either 02 or 04 in 
graph C\ to cut out one of the flags. This altered C\ then represents the untraversed remainder of C one 
must encounter, if there were no way to reach a flag without meeting both end points of a chord. □ 

Definition 6.2. An orientation-reversing q-chain Q in the simplicial complex K is an ordered chain of q- 
simplices ao, a±, . . . , o~k-i where each cxj has common p-faces tj with <7(j+i) mo( j & and tu_u mo d k with 
a (i-i) mod fc> an d the sum of the 2k entries of [d q ] indicating each n is a face of <jj and ov i+1 \ mod k is 
2 mod 4. Each Tj is an interior p-simplex of the chain. We allow g-simplices or interior p-simplices of 
Q to be repeated, as long as for any two instances of such a simplex, the simplices of the other dimension 
immediately before and after these two instances form a set of four distinct simplices. Each p-face of any 
a G Q that is not also the face of either <7u+u mo d k or a (i-i) mod fc> f° r some i indicating the order of an 
instance of a in Q, is an exterior p-simplex of Q. 

The restriction on repetition is equivalent to no entry of [d q ] being used twice in this chain. It is then 
immediate that each b-odd circuit in the bipartite graph representation of [d q ] represents an orientation- 
reversing chain, and vice versa. There is an MNTU submatrix (MNTUS) M of [d q ] whose columns are the 
columns of Q if and only if this b-odd circuit is induced, and does not properly contain another induced 
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b-odd circuit. If there is such an MNTUS, we call the rows of [d q ] that intersect M, which correspond 
to the interior p-simplices of the orientation-reversing chain, interior rows, and call the columns of [dq] 
corresponding to the g-simplices in the orientation-reversing chain Qm, and also call the rows of [d q ] that 
correspond to exterior p-simplices exterior rows. These are the rows of [d g ] that do not intersect M, but 
have nonzero entries in Qm- 

Note that if there are repeated simplices in Q, the ordering given is not unique, and two different or- 
derings may differ by more than a choice of a starting simplex Co. The repeated simplices imply that the 
bipartite graph representation is not a cycle, and a choice of ordering the simplices in Q corresponds to a 
choice of traversal of its bipartite graph. 

Definition 6.3. For a given matrix A, a columnwise minimally non totally-unimodular submatrix, or CM- 
NTUS, M of A is an MNTUS where no MNTU M' that is also a submatrix of A exists such that the set 
of columns of M' is a subset of the set of columns of M. If there is such an M', then M' is columnwise 
contained in M. 

We now describe a useful property of a CMNTUS, and illustrate the distinction between a MNTUS and a 
CMNTUS on a 2-complex in Figure 5. 

Theorem 6.4. If M is a CMNTUS of [dq], then each exterior row for M has an odd number of nonzero 
entries in Qm- 

Proof. Let % be an exterior row of an arbitrary MNTU submatrix M of [dq] with an even number of nonzero 
entries of % in Qm- Let C be the chordless b-odd circuit that is the bipartite graph representation of M. If 
C is not a cycle, split nodes as necessary to represent C as a cycle, or "wheel". Add the bipartite graph 
edges of i in Qm, and think of these edges of i as spokes of the wheel C. Call any portion of C in between 
consecutive spokes, along with these spokes, a "slice" of the wheel. Each of these slices then is a cycle, and 
there are an even number of these slices in the entire wheel. 

If C is not a cycle, then we have a choice of ordering the wheel as we split nodes to create it. To show 
our result, we choose this ordering in such a manner that there is no chord for any slice of the wheel (when 
thinking of any slice as a separate cycle). If we think of the spoke ends as being flags, Theorem 6.1 shows 
this step can be performed. 

If we start with a single slice of the wheel, and re-build the wheel by adding adjacent slices, it must be 
true that at least one of these slices must be a b-odd cycle. After putting together an odd number of b-even 
slices, the cycle that is the portion of the wheel we have built plus the two boundary spokes must be b-even. 
Hence if we put all but one of the slices of the wheel together, the resulting boundary is a b-even cycle. The 
last slice, unlike all the other previous slices, must have two edges in common with the part of the wheel 
already built. We know the entire wheel is b-odd, and so this last slice must be b-odd. We may use a similar 
b-parity argument to show that the number of b-odd slices is odd. 

Now restore C to its original form. The slice of the wheel that was a b-odd cycle is now a b-odd circuit. 
And by our choice of traversal, it is chordless. The slice contains two edges not in C, and because no two 
g-simplices may have more than one common p-face, excludes more than two edges of C. Therefore the 
submatrix M' whose nonzero entries are the edges of C contains fewer columns than M. Also, all the 
columns of M' are also columns of M. And because C is a chordless b-odd circuit, M' is an MNTU. But 
this result contradicts the assumption that M is a CMNTUS. □ 

For any MNTUS M of [d q ] with r rows, let ^# be the set of elements of Kev(A) whose nonzero q- 
coefficients are contained in Qm- Because det(M) is nonzero, Ker(M) is trivial. This means there is a 
bijection between the set of linear combinations of columns of M, and the set of possible row sums of M. 
This, along with Lemma 3.5, implies that for any m 1 , m 2 e jtft , the set of ^-coefficients of interior rows of 
m 1 and m 2 are equal if and only if all g-coefficients of m 1 and m 2 are equal. 




Figure 5: A 2-complex illustrating an MNTUS and a CMNTUS. The Mobius strip formed by all triangles 
represents an MNTUS M of the 2-boundary matrix of the complex, but M is not a CMNTUS. Edge ad 
is the face of two triangles in this Mobius strip, and thus corresponds to an exterior row that has an even 
number of nonzeros in Qm- The smaller Mobius strip formed by the seven triangles leaving out die, dhi, 
and dch represents a CMNTUS M'. Edge ad is an interior row of M' . 

Definition 6.5. For any row % of M, let [in'] denote the equivalence class of elements of ^ whose p- 
coefficients of interior rows is the unit vector with its nonzero coefficient at row i. 

Letting Pm represent the set of interior rows of M, another consequence of the kernel of M being trivial 
is the following equation for any m e ^# with p-coefficients p, and some representative m' from each [m 1 ] . 



m 



iePM 



Pirn 



(4) 



Lemma 6.6. For any MNTU submatrix M of [d q ], and any interior row i of M, each q-coefficient of any 
element of [m 1 ] is nonzero if and only if it corresponds to a column of Qm, and each such coefficient is 

±(1/2). 

Proof. For an MNTUS M of [dq], let i be an interior row. Because det(M) ^ 0, there is some entry Mij of 
M that is nonzero. If we multiply Mij by — 1, and call the result M', then the M' is Eulerian, and the sum 
of its entries is mod 4. Therefore M' is totally unimodular. Therefore there is some set ^of columns of 
M' that we may multiply these columns each by — 1, and if we call the result M'_ j , the sum of each row of 
M'_j must be 0, 1, or — 1. Because M'_ * is also Eulerian, each of these row sums must be 0. If we now 
multiply M^ by —1 again and call the result M-j-, the row sum of i is ±2, and every other row sum is still 
0. If the row sum of i is —2, multiply all columns of M_j by — 1. Now call this result, whether this final 
inversion is necessary or not, Mj. The row sum for i is 2, and all other row sums are 0, and M, is M with 
some (perhaps empty) set of columns of M scaled by — 1. Therefore by Lemma 3.5 and equation (4), any 
element of ^# whose nonzero ^-coefficients agree with the scalings of Mi is twice some m' G [m 1 ], and all 
these coefficients are either 1 or —1. Since any elements of ^# are equal in all q-coefficients if and only if 
they are equivalent, any m l G [m J ] must have all nonzero g-coefficients be ±(1/2). □ 
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From this result, the following Lemma is almost immediate, and its proof is omitted. 

Lemma 6.7. For a given MNTUS M, and any list of interior rows i\ , 12, ■ ■■i n > an d elements m' 1 , m' 2 , . ,.m'° 

o/[m !l ], [m' 2 ], ...[m ln ], respectively: 

PI. Ifn = 2, all q-coefficients of both m n + m 12 and m" — m 42 are in {0, ±1}. 

P2. Ifn = 2, for any j £ Qm> the jth q-coefficient o/m" + m 42 is zero if and only if the jth q-coefficient 
ofni 11 — m' 2 is nonzero. 

P3. Ifniseven,J2 a =i n m la is integral. 

P4. If n is odd, every nonzero q-coefficient of ^2 a= i n m J " is nonintegral, with each of these nonzero 
q-coefficients of the form | with k an odd integer 

Note that nowhere in Lemma 6.7 is it required that i a / ip for any a, /3, G 1, 2, ...n. 

Definition 6.8. For any MNTUS M of [d q ], and any concise z, let m(z) be the unique element of jtft where 
{z, m(z)} is linearly concise, Zj = implies the jth entry of m(z) < Vj < 2(m + n), and for each 
interior row i of M, the 2th p-coefficient of z and m(z) are equal. 

Theorem 6.9. For any MNTUS M of [dq], and any interior row i, there is a unique vertex z % of Pi whose 
nonzero q-coefficients are contained in Qm and whose p-coefficients of interior rows are all 0. This vertex 
is z — m(z ) where z is the identity solution to OHCPi. Furthermore, if M is a CMNTUS of [d q ], z % is 
the only nonintegral vertex of Pi whose nonzero q-coefficients are contained in Qm- 

Proof. Let z % = z 1 — m(z / ). Then z* is feasible, and all p-coefficients of interior rows of z* are zero. 
Because Ker(M) is trivial, this is the only concise feasible solution whose nonzero ^-coefficients are con- 
tained in Qm-, with the p-coefficients of interior rows all zero. To show z* is a basic solution, we try to 
decompose — m(z J ) into z c + z D satisfying the properties specified in Theorem 3.9. Because Ker(M) is 
trivial, there must be some z-coordinate that is nonzero in z l + z c , but zero in z\ implying this coordinate 
is also nonzero in z D . Therefore Property 7 in Theorem 3.9 cannot be satisfied, and so z* is a basic solution 
of OHCPj. 

Now assume M is a CMNTUS. Let z be a basic solution not equivalent to z % whose nonzero q- 
coefficients are contained in Qm- Then by Lemma 3.7, the ^-coefficient for some exterior row iq is nonzero 
in z* but zero in z. 

Since M is a CMNTUS, by Theorem 6.4, iq has an odd number of nonzero entries in Qm- Then 
Lemma 6.6 implies that an odd number of the g-coefficients for the columns with these nonzero entries must 
be integral in z. 

Decompose z into z 1 + z R + z z , where z R is the element of .^# whose nonzero (/-coefficients are equal 
to the nonintegral ^-coefficients of z, and z z is the element of ^ whose nonzero ^-coefficients are equal to 
the integral ^-coefficients of z. Then z R and z z cannot both be 0. If z 1 + z R is a basic solution, then by 
Cramer's Rule, this implies there is a non-TU matrix contained in the columns with nonzero ^-coefficients 
in z R , contradicting M being a CMNTUS. 

If z 1 + z R is not a basic solution, then by Theorem 3.3, there is some nonzero z K e Ker(^4) where 
all p- and ^-coefficients nonzero in z K are also nonzero in z 1 + z R . Since z is a basic solution, z z must 
cancel one of these coefficients. And since no nonzero ^-coefficients in z R are nonzero in z z , these canceled 
coefficients must all be p-coefficients. Since all (/-coefficients of z z are integral, all p-coefficients of z z are 
also integral. Then by Lemma 3.10, there must be at least two nonzero integral p-coefficients in z 1 + z R , 
and so there must be at least one in z R . If we construct the input chain whose nonzero coefficients are the 
values of each of the integral p-coefficients of z R multiplied by —1, then z R added to the identity solution 
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for the OHCP with this input chain must be a basic solution to this OHCP. Then by Lemma 4.3, there must 
be some OHCPj where z R added to the identity solution for OHCPj is basic. Hence, again by Cramer's 
Rule, there is a non-TU submatrix contained in the columns with nonzero (/-coefficients in z R , contradicting 
M being columnwise minimal. □ 

Lemma 6.10. For a set of columns Q of [d q ], let there be no i such that OHCPi has a nonintegral basic 
solution in (Pa); whose nonzero q-coefficients are contained in Q. Then for any basic solution z Y of any 
OHCPi in (PaJi whose nonzero q-coefficients are contained in Q, all ^-coefficients ofz are in {0, ±1}. 

Proof. If z Y is a basic solution in (Pa)^ to OHCPj for some i whose nonzero (/-coefficients are contained 
in Q, it must be integral. Let z K = z Y — z 1 , where z 1 is the identity solution. If z Y has a coefficient a 
for the ^-coordinate iq where \a\ > 1, then by Lemma 4.1 and Theorem 3.9, if we let z /o be the identity 
solution to the OHCP LP Oq whose input chain has all zeros except for coordinate Iq, which has coefficient 

±1 that is opposite in sign to the coefficient of iq in z Y , then z I(> -\ — z K is a basic solution to Oq with the 

a 
p-coefficient of i being the nonintegral value — (1/a). However, any basic solution of Oq where the nonzero 

q-coefficients are contained in Q must still be integral, giving us a contradiction. □ 

7 NTU Neutralized Complexes 

We define the concept of NTU neutralization, and present results characterizing this condition. 

Definition 7.1. For any interior row i of an MNTUS M of [d q ], let k* represent a concise integral element 
of Ker(yl) whose sum of p-coefficients of interior rows is odd, (k* — m(k 1 )) \% / 0, and let the absolute 
value of each p-coefficient of k* — m(k l ) be less than or equal to the absolute value of this coefficient in 
z\ If each interior row i of M has such a k\ then M is neutralized. If all MNTU submatrices of [d q ] are 
neutralized, then K is NTU neutralized in the qth dimension. 

Theorem 7.2. For any MNTUS M of [d q ], the projection z % \ % for each interior row i is a convex combina- 
tion ofz^\x and z 2 \%- where both z 1 and z 2 are integral elements of Pi if and only if M is neutralized. 

Proof. First, we the ignore the restriction that z 1 and z 2 are integral, and say by Corollary 5.3 that each z % \ % 
is not a basic solution to OHCPi |,sr, and hence a convex combination of z x \x and z 2 \% with z 1 , z 2 e Pi, 
if and only if there exists a k* x with k* x | % ^ 0, and where all nonzero coefficients of k* x | <% are nonzero in 
z\ 

If for a given z % there is a k* satisfying Definition 7.1, by Lemma 3.11 we may adjust k* and m(k l ) if 
necessary so that {k l , z l , m(k*)} is linearly concise. Then k* — m(k t ) satisfies the conditions for k* x . So 
let k* x = k* — m(k l ). Since k* is integral, by Lemma 6.7, each ^-coefficient in Qm of k* x is nonintegral 
with a denominator 2. By Lemma 6.6 and Theorem 6.9, this is also true of z\ Every other (/-coefficient 
in both k* x and z l is integral. So if we let z 1 = z % + k* x , and z 2 = z* — k* x , z 1 and z 2 are both integral. 
Because the absolute value of each p-coefficient of k* x is less than or equal to the absolute of this coefficient 
in z\ z 1 ! x and z 2 | x are both in Pi\%. If either z 1 or z 2 is not in Pi, then by Corollary 3.14 we may use the 
same method as in Corollary 3.13 to transform either into some solution that is concise, integral, feasible, 
equivalent, and with all x-coefficients unchanged, keeping z l \x a convex combination of z l \x and z 2 | % . 

Now suppose that for some i, z % \% is a convex combination of z 1 ) % and z 2 | % with z 1 and z 2 feasible 
and integral. Then there must be some k* x satisfying the same qualities as above, and some a > and 
(3 < where z l + ak* x and z l + /?k* x are both integral with their projections in Pi\%- 

Let a be of least absolute value such that z % + ak* x is integral. Then because both (z l + ak* x )| jr and 
(z l — ak* x )|^r are in Pj\%, the absolute value of each p-coefficient of ak* x is less than the absolute value 
of this coefficient in z\ 
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All (/-coefficients of ak* x in Qm must be nonintegral with denominator 2, and all other (/-coefficients 
integral. Hence for any interior row j, there is some m? G [m J ] where ak* x + m- 7 is concise and integral. 
Then any such ak* x + mP satisfies the conditions for k\ D 

We now present out main result, which states that the complex being NTU neutralized is equivalent to none 
of the nonintegral vertices of the OHCP LP projecting down to vertices in the projection P\%. Hence we 
cannot have a unique fractional optimal solution when the complex is NTU neutralized. 

Theorem 7.3. For a given complex K with boundary matrix [d q ], the projection z\x of each nonintegral 
vertex z of any OHCP LP over K with integral input p-chain c and polyhedron P is not a vertex of P\sc if 
and only if K is NTU neutralized in the qth dimension. 

Proof. First, we show K being neutralized is a necessary condition. If K is not neutralized, then by Theo- 
rem 7.2, for some i and MNTU M, there is a z % where z % | % is not a convex combination of z 1 1 x and z 2 1 % 
with both z 1 and z 2 integral elements of Pi . If z % | % is not a convex combination of any two points in Pi \ % , 
then it is a vertex of P\ \ sc • 

If z % | % is a convex combination of z 1 1 % and z 2 1 % , then by Corollary 5.3, and taking the notation from 
Theorem 7.2, there is some k* x and some rational number a > such that (z* ± ak* x ) \% is in Pi\ %. Since 
all variables of OHCPj are bounded, there exists a largest value of a for which this is true. Let a have this 
largest possible value. Then either (z % + ak* x ) | % or (z % — ak* x ) | % brings one of the nonzero coefficients 
of k* x | gf to 0. Suppose without loss of generality this is true for (z* + ak* x ) \%. If {z % + ak* x ) |jrisnota 
basic solution of (OHCPj)|^r, this means there is some other k* x , with some other a'. In this case, let ak* x 
represent the sum of each such ak* x , making (z* + ak* x ) | x a basic solution. 

If z % + ak* x is not integral, then by Corollary 5.12, and Lemma 5.6, there is some corresponding vertex 
z 1 of OHCPj that is nonintegral, and t}\x is a vertex of Pi\%. If z* + ak* x is integral, let ft be a nonzero 
rational value opposite in sign to a. There must be some such /? where (z* + /3k* x ) | % is in Pi\%. If there is 
a maximum absolute value for such a /?, then because z* | <% is not a convex combination of z 1 1 <% and z 2 1 %;, 
with both z 1 and z 2 integral, z* + /3k* x is nonintegral, and by a similar logic as the a case, there is some z 2 
of OHCPj that is nonintegral, and z 2 | % is a vertex of P%\x- 

If z* + ak* x is integral, and because Ker(yl) is rational, then if there is no upper bound for the absolute 
value of (3, because Ker(^4) is rational, there must be some z l + /3k* x that is integral, contradicting z*|^ 
not being a convex combination of z 1 ]^ and z 2 |^> with both z 1 and z 2 integral. Therefore K being NTU 
neutralized in the qth dimension is necessary for the projection z\%- of each nonintegral vertex z of any 
OHCP LP over K with integral input p-chain c and polyhedron P to not be a vertex of P\%. 

Now assume all MNTUS of [d q ] are neutralized. Suppose for some OHCP LP, there is a nonintegral 
vertex z where 2,\% is a vertex of P\sc- Then by Theorem 5.14, for some i, there is a nonintegral z' where 
z'|^r is a vertex of Pi\^, and Corollary 5.15 implies that there is such a z' where all nonzero ^-coefficients 
are nonintegral. We will attempt to find a minimal set Q of columns of [d q ] that contain the nonzero q- 
coefficients of such a z'. 

By Cramer's Rule, and the definition of column-wise minimal, Q must contain the columns of a CM- 
NTUS M where i is an interior row of M. But we know from Theorem 7.2 this is not enough because of 
some k* with z' ± (k* — m(k 1 )) integral. 

By Corollary 5.10, we must add q-simplices until there exist p-coefficients r and s as in Lemma 3.10. 
Hence r and s are exterior rows of M. Further, r, s, and i are interior rows of some orientation-reversing 
(/-chain C. 

By Lemma 6.10, all newly added (/-coefficients, and therefore all (/-coefficients of C, are the same in 
absolute value. So if the columns of C contain an MNTUS, the (/-coefficients of the simplices in C is a 
solution of the form z±j. Therefore Theorem 7.2 applies and we have another z l ± (k* — m(k*)) that is 
integral. Any solution basic in OHCPj must be some combination of this z % and the previous z', and so the 
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projection of any such basic solution must be a convex combination of at least two projections of these four 
integral values. 

If the columns of C do not contain an MNTUS, then by Cramer's Rule, no nonintegral vertices can be 
added. 

So we still have not found any vertex of Pi whose projection is a vertex of Pi\x- But each time we add 
g-coefficients in a minimal way to find such a vertex, we may repeat the same logic as above. Therefore no 
such vertex can exist. □ 

Remark 1 A. In either the left or right triangulation of the Example in Section 1.1 for the OHCP with input 
chain ef, z l \_%- is all edges shown in green or blue, each with a coefficient of 0.5 (see Fig. 1). m(z 7 )|x is 
the union of these green and blue edges, each with coefficient —0.5, together with edge ef with coefficient 
1. In the right triangulation, the cyan and purple chains are projections of two integral vertices z l \_%- and 
z 2 |$r, respectively. All coefficients in both of these projections are 1. Triangle adc satisfies all criteria for 
k\ Then z^ = (l/2)(zV + Ax)- Also, m(k*)k = k*| ^ - (l/2)(z 1 |, r -z 2 | x ), and z V - (k* - 

\x = z 1x- 



m(kO)Lr=z 2 | t r + (k i -m(k l ))L 



7.1 Connections to integral polytopes and TDI systems 

There exist conditions weaker than the constraint matrix being TU, which still guarantee that the LP has 
integral optimal solutions in certain cases [20, Chap. 21,22]. In particular, fc -balanced matrices define a 
hierarchy of such matrices, with TU matrices at one end [7]. The matrix A is fc-balanced for any k € Z>o 
if Aij 6 {0, ±1} and A does not contain an MNTUS with at most 2k nonzero entries in each row. If the 
constraint matrix A of the IP in Equation (1) is k -balanced, then for certain integral right-hand sides b, the 
polytope of the associated LP is integral [7]. At the same time, the poly tope of the OHCP LP is not integral 
even when the simplicial complex is NTU neutralized. Indeed, the constraint matrix A in not fc-balanced for 
any k in this case. 

A linear system A T y < f is totally dual integral (TDI) if the LP min{f T x | Ax. = c, x > 0} 
has an integral optimal solution for every c G Z m for which the minimum is finite. Every OHCP instance 
has a finite minimum, and when the complex is NTU neutralized, the OHCP LP is guaranteed to have an 
integral optimal solution. Hence the linear system A T y < f defined by the dual of the OHCP LP (3) using 
f T = [w T w T T T ] is TDI. This correspondence sheds some light on the complexity of checking if 
a given complex is NTU neutralized. The problem of checking if a linear system is TDI is coNP-complete 
[13], but could be done in polynomial time if the dimension, or equivalently, rank(^4) is fixed [9, 19]. 

8 Discussion 

Our results on MNTUS, in particular from Section 6, are specifically for such submatrices of the boundary 
matrices of simplicial complexes. NTU neutralized complexes define a class of LPs with unique structure - 
these OHCP LP polytopes may not be integral, yet, for every input chain, i.e., for every integral right-hand 
side, there exists an integral optimal solution. Our main result (Theorem 7.3) implies that when K is NTU 
neutralized, if an optimal solution of the OHCP LP is nonintegral then there must exist another integral 
optimal solution with the same total weight. If a standard LP algorithm finds the fractional optimal solution, 
we should be able to find an adjacent integral optimal solution using an approach similar to that of Giiler 
et al. [17] for the same task in the context of interior point methods for linear programming. This approach 
should run in strongly polynomial time. 

While checking whether a linear system is TDI is coNP-complete, it is not known whether a direct 
polynomial time approach could be devised to check if the simplicial complex is NTU neutralized. Another 
interesting question is whether the definition of the complex being NTU neutralized could be simplified for 
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low dimensional cases, which could also be tested efficiently. On a similar note, are there special classes 
of 2-complexes that are guaranteed to be NTU neutralized? Two potential candidates for such classes are 
2-complexes embedded in M 3 whose first homology group is trivial, and 2-complexes embedded in M 3 in 
which each edge is the face of at least 2 triangles. The NTU neutralized complex (in right) in Figure 1 
illustrates a case where the same orientation-reversing (/-chain neutralizes all relevant elementary chains. A 
characterization of the structure of such complexes could also prove very useful. 
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